Exchange anisotropy, disorder and frustration in diluted, predominantly 
ferromagnetic, Heisenberg spin systems 
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Motivated by the recent suggestion of anisotropic effective exchange interactions between Mn 
spins in Gai-^Mn^ As (arising as a result of spin-orbit coupling), we study their effects in diluted 
Heisenberg spin systems. We perform Monte Carlo simulations on several phenomenological model 
spin Hamiltonians, and investigate the extent to which frustration induced by anisotropic exchanges 
can reduce the low temperature magnetization in these models and the interplay of this effect with 
disorder in the exchange. In a model with low coordination number and purely ferromagnetic 
(FM) exchanges, we find that the low temperature magnetization is gradually reduced as exchange 
anisotropy is turned on. However, as the connectivity of the model is increased, the effect of 
small-to-moderate anisotropy is suppressed, and the magnetization regains its maximum saturation 
value at low temperatures unless the distribution of exchanges is very wide. To obtain significant 
suppression of the low temperature magnetization in a model with high connectivity, as is found 
for long-range interactions, we find it necessary to have both ferromagnetic and antiferromagnetic 
(AFM) exchanges (e.g. as in the RKKY interaction). This implies that disorder in the sign of the 
exchange interaction is much more effective in suppressing magnetization at low temperatures than 
exchange anisotropy. 

PACS numbers: PACS: 75.10.Nr, 75.10.Hk, 02.70.Tt, 75.60.-d 



I. INTRODUCTION 



It has recently been suggested that frustration ef- 
fects may be important for the magnetic properties 
of diluted magnetic semiconductors (DMSs) such as 
Gai^aMn^Asi 1 : 2 : 3 The reported Curie temperatures (T c ) 
in these compounds continue to rise, with a maxi- 
mum of 140K recently reported for Gai-^Mn^As^ with 
even higher temperatures reported for related materials. 5 
Aside from increasing T c , it will also be important to 
have a thorough understanding of the different aspects of 
their magnetic properties over a wide temperature range 
in order to be able to construct optimized spintronic de- 
vices. In theoretical analyses, Zarand and Jankoi showed 
that within the RKKY approximation, a proper treat- 
ment of the spin-orbit coupling leads to anisotropic ex- 
changes between Mn spins. Using this interaction they 
found that the saturation magnetization is reduced by 
up to 50 % at low temperatures. Experimentally it has 
been observed that in many DMS, the saturation mag- 
netization at low temperatures is not as large as would 
be expected if all Mn moments were aligned (i.e. the full 
saturation value) & While it is likely that magnetically 
inactive Mn, such as Mn interstitialsji may account for 
some of the suppression of the low-temperature magne- 
tization, the results of Zarand and Janko suggest that 
anisotropic spin interactions may also play a significant 
role in accounting for it. We carefully investigate this 
possibility in this paper. 

In a system consisting of localized spins (local mo- 



ments) coupled to non- interacting fermions (carriers), 
where the spin-carrier interaction is a perturbation on 
the fermion Hamiltonian, the low-lying spin excitations 
can be described in terms of an RKKY interaction be- 
tween the spins. In the case of DMS, the carrier density 
is low (in fact, lower than the density of local moments by 
a considerable factor, due to carrier compensation), and 
the carrier-spin interaction is quite strong (necessary to 
enable a high T c ). Consequently, the RKKY approxima- 
tion appears unlikely to be appropriate for a quantitative 
description of the ferromagnetism in Gai-^Mn^As since 
the Fermi energy is not necessarily much larger than the 
magnetic coupling. However, the qualitative prediction 
of anisotropic exchange may be present in more precise 
treatments £ 

In what follows, therefore, we assume that the mag- 
netic properties of the system can be modeled in terms 
of an effective spin Hamiltonian, with a form that main- 
tains the symmetry properties of the RKKY interactions 
obtained in the weak-coupling limit. Rather than cal- 
culate the anisotropy in the Mn-Mn effective exchange 
within a microscopic model, we use a different approach 
to investigate the relevance of anisotropy for systems 
of diluted spins. We consider phenomenological mod- 
els with disordered, anisotropic exchanges between clas- 
sical Heisenberg spins placed randomly at low densities 
on a fee lattice (corresponding to the Ga fee sublattice 
in (Ga,Mn)As), and study the magnetic properties for 
several functional forms of the exchange interactions. In 
each case we consider various values of the disorder and 
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anisotropy. We focus on whether full saturation in mag- 
netization is reached at low temperatures and also in- 
vestigate the magnetic susceptibility, since this is known 
to be a good experimental indicator of spin freezingiiS 
Our results allow us to infer parameter ranges in which 
anisotropy and/or disorder are likely to play an impor- 
tant role in the magnetic properties of such systems. 

The paper— is organized as follows: in Sec. |n] we 
describe the models we study and how disorder and 
anisotropy are incorporated into each of them. Sec- 
tion ^0] lists the quantities we calculate and how we per- 
form the Monte Carlo simulations. Section llVl shows our 
results for the magnetization, susceptibility and Curie 
temperature and finally in Sec. we discuss our results 
and their implications for modeling III-V DMS. 



II. MODEL 

We consider A^ spins randomly distributed at loca- 
tions Rj , i = 1 , . . • , Nd on a fee lattice of size N x 
N x N, corresponding to an impurity concentration 
x = Nd/4:N 3 . The spins are treated as classical vari- 
ables (Mn spins in Gai-^Mn^As have S = 5/2, so this 
is a reasonable approximation). For simplicity, we take 
the classical spins to have unit length; this is equivalent 
to rescaling the exchange from J to JS 2 , i.e. changing 
the units of energy. The most general formulation of 
the problem we consider in this study is provided by the 
Heisenberg Hamiltonian 







(1) 



i,j a/3 



where a and (3 index Cartesian coordinates. The ex- 
change coupling JijF a p is written as a product of a ran- 
dom variable Jij , and a function F of the separation of 
the two spins, for reasons that will become clear later. 
Since the sites are not on a regular lattice, the sum- 
mation over site index is always over all sites of the 
system. For uncoupled spins i and j we simply have 
F a p(Ri — Rj) = 0. 

We consider the exchange integral to be parameterized 
as: 



F a /3(Ri - Rj) = ( \5 aP + (1 - A)eg.e£. 



(2) 



|R; — Rj| and the unit vector e.^ = (R; 



where r 

Rj)/|Ri — Rj|. The parameter A controls the exchange 
anisotropy, where e.y defines the axis of anisotropy 
for the pair of spins Si and located at positions 
Hi and Rj. Using Eq. J3J), the effective interaction 
can now be rewritten as JijSfFapCRi — Hj)Sj — 

•WW ( S iSj +AS' i L S , i L ), where the parallel and per- 
pendicular components are defined with respect to ey. 
As a result, for A = 1 the model has no anisotropy and 
reduces to a simple disordered Heisenberg model. For 
< A < 1 the couplings become anisotropic and favor 



alignment of each pair of spins along their positional axis 
ey, whilst for A > 1, the couplings favor alignment per- 
pendicular to this axis (see Fig. Since the relevant 
directions ey differ considerably from pair to pair, some 
frustration is introduced into the system, possibly leading 
to spin-glass physics. 



^ 1 Interactions favour * 
alignment perpendicular 



to line joining spins 



X> 1 



X< 1 



Interactions favour alignment along 
line joining spins 

FIG. 1: Schematic representation of the two different types 
of anisotropy parametrized by A. If A > 1, the anisotropy 
defines an easy plane for the interactions between two spins, 
whilst if A < 1 the anisotropy defines an easy axis. 

The function /(r) describes the spatial variation of the 
exchange interactions. In this work we investigate the 
following possibilities: 



9(R C - r) Short range FM, model A,±I 



sin y — y cos y 
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RKKY, model B, 



Long range FM, model C, 



(3) 



where y = 2kpr, and Uf is the Fermi wavevector. In 
the short range model, equal-strength ferromagnetic in- 
teractions exist only between neighboring spins within 
distance R c of one another. 11 The RKKY model allows 
interactions between all spins in the system, with oscil- 
lating sign depending on the Fermi wave-vector kp cor- 
responding to different carrier concentrations. Finally 
we consider a purely ferromagnetic model with long- 
range (power-law) interaction. A comparison between 
the short-range and the long-range models clarifies the 
role played by the range of the exchange, while a com- 
parison between the RKKY and the long-range model 
clarifies the importance of the exchange sign oscillations. 

In the RKKY approach of Zarand an Jankoji the rel- 
ative magnitudes of the exchanges parallel and perpen- 
dicular to the line joining two Mn spins, K par {r) and 
K perp (r) respectively, depend on the distance r between 
the two spins. As a rough guide to compare with the 
models we use in this work, \K perp (r)\ > \K par (r)\ for 
large r while \K par {r)\ > \K perp {r)\ for small r— It should 
be noted that we use the RKKY interaction for illustra- 
tive purposes only, since the exchange interactions in a 
more realistic model of DMS will likely be quantitatively 
different. 
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As demonstrated in mean-field and Monte Carlo stud- 
ies of an impurity band model for III-V diluted magnetic 
semiconductors ^i 12 ' 1 ?i 14 ' 1 ^ 1 ^' 1 Ti 1 ?i 1 ? and in studies of the 
kinetic-exchange model of III-V DMS including disorder 
and Coulomb interactions of the charge carrier with the 
Mn acceptors^ inhomogeneity induced by positional dis- 
order of the Mn spins implies different local carrier charge 
densities at different sites, which in turn leads to a broad 
distribution of effective local fields created at various Mn 
sites by the itinerant carriers. If one integrates out the 
fermionic degrees of freedom and formulates the problem 
in terms of effective exchanges between the Mn spins, 
there should also be a wide distribution of their effec- 
tive exchanges. To incorporate the effect of positional 
disorder leading to a wide distribution of magnitudes of 
the exchanges in our model, we assume that Jij — ^z^j ; 
where the are random variables attached to each site, 
such that their logarithm Zi = log 10 e-j has a Gaussian 
distribution 



"PnM) = 



1 



27TCT 



exp 



(zj - z{ni)) 2 
2a 2 



(4) 



although the precise form of the disorder should not 
change our conclusions. In Eq. Q, the mean value of 
Zi, Zi(rii) is taken to depend on the number of nearest 
neighbors spins (i.e. within a distance R c ) of the site i, 
Hi though the relation z~i(rii) = a(rii — no), where no is 
the average number of neighbors for a given cutoff R c . 
In our simulations R c is chosen such that no — 6 or 12. 

This scheme naturally favors stronger couplings within 
more dense clusters (higher rij values) and hence should 
mimic some of the phenomenology observed in the impu- 
rity band model for III-V DMS, where holes congregate 
in regions of higher Mn density and therefore lead to 
stronger effective interactions between the Mn in these 
regions^ The parameter a controls the width of the dis- 
tribution of couplings, and therefore we use it to charac- 
terize the disorder present in the system; a more disor- 
dered system with a wider distribution of couplings cor- 
responds to a larger value of a. 

Thus, the three parameters that control the behavior 
of the model are: (i) A, which controls the amount of 
anisotropy; (ii) cr, which controls the disorder-induced 
width of the distribution of effective exchanges, and (iii) 
R c , which defines the average number of nearest neigh- 
bors for the short range model (for the RKKY and long 
range models, R c = oo). 



A. Parameters and calculated quantities 

We investigated systems of linear sizes L = 11, 14 and 
17, containing Nd — 53, 110 and 196 spins respectively. 
This corresponds to a Mn density in Gai-^Mn^As of 
x = 0.01. The interaction range R c is chosen such that 
the average number of nearest neighbors within R c is 
either no = 6, or no = 12. We studied the short range 
model for both values of no- Note that for models with 
short-range interactions, the value of x has no qualitative 
effect on the behavior of magnetic properties for x < 
0.05, since changing the concentration is equivalent to a 
pure rescaling of all inter-spin distances by a fixed factor. 
The value of x becomes relevant for larger concentrations, 
where the average inter-spin distance is comparable with 
the lattice constant. The anisotropy values considered 
were in the range A = 0.1toA = 10, with the isotropic 
case A = 1 used as a reference. 

For the short range model, two different disorder val- 
ues were considered for most anisotropy values, whilst for 
the other models, only one disorder strength was consid- 
ered. The values of a considered were: in the short range 
model, a = 0.017 and 0.05 for no = 6, and a = 0.01 and 
0.03 for no = 12. In the long range model we chose 
a = 0.03, while for the RKKY model, a = 0.01. These 
values were chosen such that meaningful comparisons be- 
tween different models can be performed, subject to some 
computational constraints. Figure|2]shows the actual dis- 
tribution of Jij used in the simulation for two different 
values of a. 




III. SIMULATIONS 

We have performed exhaustive Monte Carlo simula- 
tions on each of models A, B and C for temperatures 
above the ordering temperature T c to well below T c . We 
use a range of sizes to determine T c via finite size scal- 
ing, and average over a sufficient number of realizations 
of the disordered systems for each size to determine av- 
erage quantities. 



FIG. 2: Distribution of log 10 Jij. Each curve is an average 
of 5 samples with L = 32, Nd = 1310 and no = 6. The 
distribution is almost Gaussian, and its width increases with 
increasing a. 

For each pair of A and a values we considered many 
realizations of positional disorder of the Mn spins (gen- 
erally at least 40) and then averaged over these configu- 
rations to obtain our final results. 

We calculated equilibrium (disorder averaged) aver- 
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ages for the following quantities: (i) the magnetization 



C. Monte Carlo Technique 



M 



5> 



(5) 



and (ii) the magnetic linear susceptibility 

2v 



Xrn = (3 




M 



f2 



(6) 



To determine the Curie temperature in our samples we 
calculated the Binder cumulant 



/ 



G(L,T) = - 



5-3- 



V 



~m Si s > 



(7) 



and used finite size scaling. G(L, T) is defined such that 
in the paramagnetic phase it decreases with L, and tends 
to zero as L — > oo, while in the ferromagnetic phase it 
increases with increasing size L and tends to unity in 
the thermodynamic limit. Near the transition temper- 
ature T c , this dimensionless quantity has the finite size 
scaling form G(L,T) = G[L~ (T - T c )\ where v is the 
exponent of the diverging spin-spin correlation length 
f ~ (T - Te)^, 21 ! 22 ! 23 Consequently, at T c , G(L,T C ) is 
independent of L; T c can be identified by a simultaneous 
crossing of G(L, T) vs. T curves for different L. This 
method is found to be more reliable in determining T c 
than the onset of magnetization, or the position of peaks 
in the magnetic susceptibility in relatively small finite 
size samples^. 



The Metropolis algorithm necessitates long equilibra- 
tion times at low temperatures and large values of 
anisotropy, due to its relatively slow sampling of phase 
space. As a result, we used a different method to per- 
form spin flips. This method is formally equivalent to the 
Metropolis algorithm and can be used for any Hamilto- 
nian, TL, involving classical spins that can be written in 
the form Tt = hi • Si, where h; = Y2j JijSj is the local 
field created at site i by the other spins. In our model, 
it also contains terms of the form y\ - Jjjeij(ey • Sj). 

The implementation of a Monte Carlo simulation re- 
quires successive spin flips at each site in the system for 
a Monte Carlo step. Each spin flip involves changing the 
angular position vector of one spin, whilst keeping all oth- 
ers fixed. After a sufficiently large number of spin flips at 
each site, the angular distribution of these vectors will be 
equal to the equilibrium (Boltzmann) distribution. For 
any spin that is about to be flipped, the distribution of 
the angle 9 between the spin and its local field (inte- 
grated over azimuthal angle) is 



P{0) = 



ke 



suit 



,-k 



(8) 



where k = /3|hj|. This distribution is easily integrated 
over 9, and hence by finding a function f(x) which maps 
x £ [0,1] to cost? £ [—1,1], uniform sampling of the 
[0, 1] interval will give the desired distribution of cosQ £ 
[—1, 1]. It is simple to find that f(x) is given by 



,9 

f- l (cos6) = x= / p{6')d6', 
Jo 



and the explicit expression for this mapping is: 



(9) 



B. Temperature Rescaling 

In order to compare temperature scales for differ- 
ent A values, we note that the exchange felt for a 
spin orientation at angle 9 to the axis is J{9) = 

J V cos 2 9 + X 2 sin 2 (9. Thus, (J) = ± j dttJ(9) gives 
the average exchange integrated over all solid angles fi. 
Evaluating the integral gives^ 



(J) = Jo x 



1 



1 

i sinh~ 



i ( VEE£k 



A 



,A> 1, 
,A = 1, 
,A < 1. 



These factors are used to rescale the temperature for 
each value of A chosen. The plots are for the temper- 
ature scaled as T/J c g(X), where J e ff(l) = ^cfr(O-l) = 
0.5150 J , J e ff(0.5) = 0.6901 J , J e ff(1.5) = 1.3463 J and 
J eff (10) = 7.8902 J . 



cos6 = f(x) = l + -hi[l-x(l-e- 2fc )]. (10) 
k 

This approach is similar to those that are sometimes used 
in lattice field theories. 25 We have tested that our method 
yields results for equilibrium quantities that are identi- 
cal to the Metropolis algorithm, at the same time, the 
equilibration and auto-correlation times were found to 
be shorter by a factor from 10 to 100 depending on tem- 
perature and size when using this method for the above 
models for our range of T and L. 

If the new orientation of the spin 9 is generated in 
this way, each trial move is accepted, since the mapping 
we have made is one that changes the probability of a 
configuration being accepted from the Boltzmann weight 
to unity, provided the configuration is chosen in the pro- 
scribed manner. The calculation of the mapping function 
may be more time-consuming than calculation of the en- 
ergy difference between initial and final states performed 
in the Metropolis algorithm, but it has the advantage of a 
100% acceptance rate. This compares to the Metropolis 
algorithm which may have exponentially low acceptance 
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10 100 1000 

Monte Carlo Steps 



FIG. 3: A comparison between Metropolis and our algorithm 
shows the convergence to equilibrium of Monte Carlo sim- 
ulation is accelerated about 10 times. This simulation was 
done with 196 spins, T = 0.5, A = 0.5 and a = 0.1. Each 
MC step flips 196 spins randomly. Equilibrium is deter- 
mined by convergence between replicas started with all spins 
aligned(ferromagnetic initial configuration), and random ini- 
tial configuration. At lower temperatures and with larger 
disorder, an increase in speed of more than a factor of 100 
was observed. 



rates. The evaluation of the mapping function can be 
optimized to further to increase the efficiency. The cal- 
culation of k takes about the same amount of time as 
evaluating the energy difference in the usual implemen- 
tation of the Metropolis algorithm. After updating cos0 
all that remains is to generate a random azimuthal ori- 
entation of the spin by randomly selecting a vector (of 
appropriate magnitude) perpendicular to the local field; 
this can be done with a fast algorithm. 

In summary, the algorithm used is: 1) compute the 
local field and k = /3|hj|; 2) choose a random number 
x £ [0, 1] and use Eq. I|1U|) to generate the new orienta- 
tion cos#; 3) update the energy; 4) generate a random 
azimuthal component to the spin; 5) update spin compo- 
nents; 6) repeat for the next spin. 

We have compared this algorithm with the Metropolis 
algorithm for our models. They generate the same equi- 
librium results, however this algorithm equilibrates much 
faster when the model has large anisotropy and frustra- 
tion. Figure |3 shows an example in which the algorithm 
we used reached equilibrium an order of magnitude faster 
than the Metropolis algorithm. 



IV. RESULTS 

We present the results we obtained for the magnetiza- 
tion for all three models. For the short range model, we 
also present the susceptibility and the Curie temperature. 



A. Short range model 

1. Magnetization 

The magnetization curves M (T) shown in Fig. 0] for 
size corresponding to Nd — 196 spins have the character- 
istic linear decrease with temperature seen in previous 
work and in experiments in Gai-^Mn^Asi^SiSi It is also 
apparent from comparing Fig. EJa) with Fig. 0Jb) that 
increasing the number of nearest neighbors in the short 
range model does not qualitatively change the shape of 
the magnetization curve. 

With temperatures properly rescaled by J c ff, we ob- 
served that curves corresponding to the same disorder 
value a but various anisotropies A are almost identical at 
high temperatures. Increasing the disorder a leads to in- 
creased magnetization at high temperatures, as expected 
from studies of other models for DMSii^iii^ 

Anisotropy plays a role at low temperatures, where in- 
creased anisotropy does lower the T = magnetization, 
although the effect is rather small. The suppression is 
more pronounced for A < 1 (easy-axis) than A > 1 (easy- 
plane). For no — 6, the value of the saturation magne- 
tization seems to depend only on A and be independent 
of a. Whereas at higher temperatures (T > 0.5 T c ), as 
stated above the opposite is true. Consequently, both 
anisotropy and disorder in the magnitude of the fer- 
romagnetic exchanges affect M(T) in the temperature 
range T = 0.1 T c ~ 0.5 T c . For n Q = 12, the suppres- 
sion of low-temperature magnetization by anisotropy is 
not observed except for an extremely large anisotropy 
A = 0.1. 



2. Susceptibility 

The behavior of the linear susceptibility for the short 
range model with both 6 and 12 nearest neighbors is 
shown in Fig. [3] After temperature rescaling, the high- 
temperature tails of the curves with the same a value are 
again identical. At low temperatures, anisotropy leads 
to a finite value for the T = linear susceptibility. The 
finite value of the linear susceptibility at T — is consis- 
tent with the incomplete saturation of magnetization in 
the presence of anisotropy. This effect is most transpar- 
ent in the case of no = 6; in contrast, for tiq = 12, the 
effect of anisotropy is very weak unless the anisotropy is 
extreme, e.g. A > 5 or A < 0.1. 



3. Curie Temperature 

The Curie temperatures are deduced from the Binder 
cumulant Eq. @ curves. Figure shows curves cal- 
culated for different system sizes from the 12 neighbor 
model which intersect at the same temperature for a 
given A, indicating the position of T c . These critical 
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FIG. 4: Magnetization for (a) no — 6 and (b) no = 12, for 
various values of A and a. The simulations are for Nd = 196 
spins. The insets amplify the low temperature regions. 



FIG. 5: Magnetic susceptibility for (a) no = 12, and (b) 
no = 6, with various values of A and a. The simulations are 
for Nd = 196 spins. 



temperatures are very close to the peaks in linear sus- 
ceptibility of the same model, and are tabulated in Table 
H] Note that the temperature is not rescaled by J c s in 
Fig. 0in order to give a clear view of each set of curves. 
When the temperature is rescaled, the values of T c are 
approximately the same. 



Tc/J e ff 


A = 0.5 


A = 1 


A = 1.5 


From xm 


3.42 


3.61 


3.50 


From G(N, T) 


3.39 


3.48 


3.39 



TABLE I: Critical temperatures (in units of Jeff, see text) 
estimated from Binder cumulant G(N, T) and magnetic sus- 
ceptibility \M f° r a — 0.01 and no = 12. 

We conclude that small-to-moderate anisotropy has 
a rather small effect on the magnetic properties of 
the short-range model. Significant deviations from the 



isotropic behavior are seen only for very large values of 
the anisotropy, and are generally more pronounced in the 
no = 6 model. 

B. Long range model 

The qualitative behavior of the low-temperature mag- 
netization in the long range model is very similar to that 
observed in the short range model, as can be seen in 
Fig. [7| Increasing the length-scale of the exchange does 
very little to increase the amount of frustration in the 
model, and in both the short-range and the low-range 
model the suppression of the low temperature magneti- 
zation is less than 20%, and often much smaller, for all 
values of anisotropy considered. This is small compared 
to the 50-60% reduction observed experimentally, 6 hence 
other sources of frustration must be present to account 
for the greatly reduced low temperature magnetization as 
compared to that of an aligned ferromagnet. One candi- 
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FIG. 6: Binder cumulant for models with 12 neighbors on 
average, a = 0.01, various anisotropy A and sizes L. 



FIG. 8: Magnetization for the isotropic, A = 1 RKKY model, 
at various values of the Mn concentration x and hole compen- 
sation 1 — p (corresponding to different choices of y = 2/cft). 



date for the source of further reduction is the antiferro- 
magnetic component of effective Mn-Mn interactions as 
is found for RKKY interactions. 
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FIG. 7: Magnetization for no = 6, various A, a = 0.03, long 
range model. In this case all pairs of spins interact, and no is 
used only to adjust the average coupling constant. 



C. RKKY models 

Figure [S] shows the magnetization calculated for 
isotropic RKKY models with different Fermi wave- 
lengths, corresponding to different charge carrier con- 
centrations. The main observation is that the antifer- 
romagnetic exchange at long distances introduced by the 
RKKY model has a significant effect on lowering the low 
temperature magnetization; the decrease is considerably 
larger than that due to anisotropy in purely ferromag- 



netic exchange interactions in the short- and long-range 
models considered previously. 

In Fig. IS] we show curves corresponding to several dif- 
ferent values of the Mn concentration x and charge car- 
rier compensation 1 — p; the resulting charge carrier con- 
centration is n/i = 4xp/a 3 , where a — 5.65 A is the 
GaAs lattice constant. Note that the temperature also 
needs to be rescaled by a factor of (px) 3 in 3 dimensions 
for RKKY interactions^ to compare different magneti- 
zation curves. One interesting observation of the series 
with x — 0.05 is that increasing the hole concentration 
px at with fixed x does not appear to change T c but 
the low temperature magnetization has a minimum at 
about p = 0.3. Further increments of p increase the low 
temperature magnetization gradually. There is a large 
range of values of the T = magnetization that can 
be achieved by tuning the Mn and the hole concentra- 
tions. This can be understood qualitatively in the fol- 
lowing way: larger hole concentrations px lead to larger 
Fermi wavevectors hp. As a result, the oscillation be- 
tween ferromagnetic and antiferromagnetic interactions 
described by the RKKY interaction appears at a shorter 
distance. If this distance becomes comparable to the av- 
erage inter-spin distance, significant frustration is present 
in the system leading to suppressed low-T magnetization. 
When p = 0.3, the antiferromagnetic part seems to be 
dominant. The magnetization actually decreases as tem- 
perature is lowered below 0.4 J e g. This effect arises due 
to antiferromagnetic coupling between clusters as a result 
of the long-range nature of the RKKY interaction^ 

We have also investigated RKKY models with ex- 
change anisotropy as discussed before. In all cases, the 
additional suppression of the low-T magnetization in- 
duced by the anisotropy is very small (less than 10%). 
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V. DISCUSSION 

The major result that can be deduced from our simu- 
lation is that while exchange anisotropy can change the 
temperature scale at which ferromagnetism occurs in a 
model, it does not change the essential qualitative fea- 
tures such as the shape of the magnetization curve, or 
the form of the linear or susceptibility greatly above T c , 
unless the ratio between the parallel and perpendicu- 
lar exchanges is extremely different from 1. Below T c , 
the anisotropy preserves a finite magnetization and lin- 
ear susceptibility at T = 0. This effect is suppressed by 
higher connectivity (large no). [Locally the random lat- 
tice with no = 12 is more symmetric than for no = 6, so 
that the anisotropy a spin feels from one neighbor is more 
likely to be cancelled by the effects of other neighbors.] 
However, in general the effect of anisotropy is relatively 
weak, as it does not lower magnetization at T — by 
more than 20%. The presence of antiferromagnetic in- 
teractions, as occur in oscillatory exchange interactions, 
such as in the RKKY model, appears to be much more 
important for lowering the magnetization at low temper- 
ature. Thus the large reduction in the magnetization 
seen in Ref. [l| can probably be attributed to RKKY in- 
teractions, rather than the effects of anisotropy in the 
exchanges. 

Our result on the effects of moderately anisotropic ex- 
change interactions on the magnetic properties of DMS 
are in substantial contrast with the results of Zarand 
and Jankoi They found large differences between the 
low temperature magnetization (but similar Curie tem- 
peratures) in the isotropic and anisotropic cases. We are 
puzzled by this difference, but suggest that it may be due 
to the fact that they used an exponential cut-off in their 
exchange interaction. This cut-off may have damped the 
effect of anti-ferromagnetic interactions in the isotropic 
case, but since the parallel and perpendicular exchanges 
had different spatial dependences in their study, they 
may not have reduced the effect of anti-ferromagnetic 
exchanges in the anisotropic case. 

There are a number of other interesting observations. 
The two types of anisotropy, corresponding to the cases 
A < 1 and A > 1 are qualitatively different. For A < 1, 
the spins preferentially align along the line joining them, 
whilst for A > 1 there is a preferred plane in which the 
spins may lie. This difference in dimensionality appears 
to explain why the magnetization is not suppressed as 
much at low temperatures for A > 1 as compared to 
A < 1, since it is easier for spins to relax from frustration 
in two rather than one dimensions. 

Whilst anisotropic exchange interactions are one pos- 
sibility which can lead to a reduction in the saturation 
magnetization at low temperatures, there are a number of 
other possibilities. Firstly, there are believed to be anti- 
ferromagnetic nearest neighbor interactions between Mn 
spins, that should have little effect due to the diluteness 



of the Mn spins, but could lead to a decrease in the sat- 
uration magnetization. 3!l Next is a theoretical proposal 
that there is an instability purely in the presence of dis- 
order towards a non-collinear ground stateASi Another 
is the experimental observation that there are significant 
numbers of interstitial Mn that appear to be involved in 
compensation processes and thus do not polarize. 7 It is 
not completely clear to what extent the saturation at low 
temperatures in Gai-^Mn^ As is due to non-participation 
of Mn spins due to the presence of Mn interstitial de- 
fects. However, given the recent progress in this area^ it 
will probably not be long before an accurate estimate of 
the proportion of Mn that are participating in the ferro- 
magnetism is known. When this is quantified, it should 
be possible to determine whether anisotropy need be in- 
cluded in realistic models of DMS and if so, how much. 
The carrier-mediated nature of the ferromagnetism also 
suggests that anisotropy may be important for transport, 
especially in the insulating phase, since hopping between 
sites will be preferred when the Mn spins have similar 
orientations i22 

We mention in passing that while anisotropy has 
small effect on the magnetization, it seems to have sig- 
nificant effect on the non-linear susceptibility at low 
temperature^ This suggests that experimental investi- 
gations of the non-linear susceptibility might shed light 
on the magnetic state in DMS. It was recently suggested 
that for x < 0.01 and x > 0.1 in Gai-^Mn^As, there 
may be a spin glass phased Here we suggest that if there 
are strongly anisotropic exchange interactions, signatures 
might be seen in co-existence with ferromagnetism. If 
such spin-glass signatures are not seen in (Ga,Mn)As, 
they may be present in other insulating materials with 
lower carrier concentrations such as Ge:Mni2i 

In conclusion, there are still a number of outstanding 
questions as to the nature of the ferromagnetic state in 
DMS, and whether anisotropy plays an important role 
in these materials. This work should be of help in clar- 
ifying which types of models are likely to be affected by 
anisotropy and what types of experimental probes might 
help to detect it. 
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